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ABSTRACT 

Based on a cosmological N-body simulation we analyze spatial and kinematic alignments of satellite 
halos within six times the virial radius of group size host halos (rvir). We measure three different 
types of spatial alignment: halo alignment between the orientation of the group central substructure 
(CCS) and the distribution of its satellites, radial alignment between the orientation of a satellite and 
the direction towards its GCS, and direct alignment between the orientation of the GCS and that of 
its satellites. In analogy we use the directions of satellite velocities and probe three further types of 
alignment: the radial velocity alignment between the satellite velocity and connecting line between 
satellite and GCS, the halo velocity alignment between the orientation of the GCS and satellite 
velocities and the auto velocity alignment between the satellites orientations and their velocities. We 
find that satellites are preferentially located along the major axis of the GCS within at least 6 rvir 
(the range probed here). Furthermore, satellites preferentially point towards the GCS. The most 
pronounced signal is detected on small scales but a detectable signal extends out to ^ 6 rvir. The 
direct alignment signal is weaker, however a systematic trend is visible at distances < 2 rvir. AH 
velocity alignments are highly significant on small scales. The halo velocity alignment is constant 
within 2 rvir and declines rapidly beyond. The halo and the auto velocity alignments are maximal 
at small scales and disappear beyond 1 and 1.5 rvir respectively. Our results suggest that the halo 
alignment reflects the filamentary large scale structure which extends far beyond the virial radii of the 
groups. In contrast, the main contribution to the radial alignment arises from the adjustment of the 
satellite orientations in the group tidal field. The projected data reveal good agreement with recent 
results derived from large galaxy surveys. 

Subject headings: dark matter — galaxies: clusters: general — galaxies: kinematics and dynamics — 
methods: numerical 



1. INTRODUCTION 

Over the last decades observational and numerical 
evidence has substantiated the picture of a filamentary 
large-scale structure in the universe. In principle the 
large-scale tidal field is expected to induce large-scale 
correlations between the orientations of halos and 
galaxies th at are embedded within these filament s (e.g.. 
F'en et al.l 2 000: Croft & Metzler 2000; Heave ns et all 
20001: ICatelan et al.l 120011: ICrittenden et all l2001t 
Porciani et al.ll2002t lJingll2002f ). On the other hand, the 



subsequent accretion onto larger systems, such as groups 
and clusters of galaxies, may alter the orientations of 
these (sub-)struc tures in response to t he local tidal field 
(|Ciotti &: Duttalll994l:lLee et al.|[2005f) . Cosmological N- 
body simulations provide a valuable tool to differentiate 
the various contributions to the halo/galaxy alignments 
within overdense regions. 

Observationally, various types of alignment between 
galaxies and their environment have been detected on 
a wide range in scales, from super-cluster systems down 
to the distribution of the satellite galaxies in our Milky 
Way (MW). On cluster scales various different types of 
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alignment are discussed in t he literature : ahgnment be 
tween neigh boring cluster s (lBing gelini982!: 'Ulmer et al 
IT989t [1981 iPlionid |199i iChambers et al. 2002^ 

between brightest cluster galaxies (BCGs) and their 
parent clusters (Carter & Metcalfe 1980; Binggeli 198_ 
Strubld [l990i : fPhee fc Latouijfl991l: IPhonis fc Basilakoi 



2002|) . between the orientatio n of satellite galaxies and 



the o rientation of the cluster (jDekell 119851 : iPlionis et al.l 
l2003l ). and between the orientation of satellite galaxies 
and the orientation of the BCG (Strublc 1990). Accord- 
ing to these studies the typical scales over which clusters 
reveal signs for alignment range up to 10 — 50 /i~^Mpc, 
which can be most naturally explained by the presence 
of filaments. 

With large galaxy redshift surveys, such as the 
two-degree F ield Galaxy Redshift Survey (2dFGRS, 
IColless et all |2001|) and the Sloan Digital Sky Survey 
(SDSS, l^rk et al.l [2000l ). it has recently also become 
possible to investigate alignments on group scales using 
large and homogeneous samples. This has resulted in ro - 
bust detect i ons o f var ious alignmen ts: Braincrd (2005), 
lYang et all (|2006( l and lAzzaro et al.l (2007) all found that 
satellite galaxies are preferentially distr ibuted along the 
major axes of their host galaxies, while [ Pereira fc KuhnI 
(2005) and I Agustsson fc BrainerdI ()2006af ) noticed that 
satellite galaxies tend to be oriented towards the galaxy 
at the center of the halo. 

In contradiction to the studies above, iHolmber^ 
(|1969| ) found that satellites around isolated late type 
galaxies preferentially lie along the minor axis of the 
disc. Subsequent studies, however, were unable to con- 
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firm th is so-called ' Holmb er g effect '(iHawlev fc PeeblesI halos. 



1971 iSharpetalJ 



t lolmb erg 



MacGiinvrav et alJ Il98l 



Zaritskv et al.l 119971 ). Recently lAgustsson &: BrainerdI 
( 2007f ) reported a Holmberg effect at large projected 



distances around blue host galaxies, while on smaller 
scales the satellites were found to be aligned with the 
major axis of their host galaxy and iBailin fc Steinmet3 
(|2005l ) claim that a careful selection of isolated late-type 
galaxies reveals the the tendency for the satellites to 
align with the minor axis of the galactic disc. I nvesti - 
gating the companions of M31 iKoch fc Grebe] ()2006D 
find little evidence for a Holmberg effect. Yet, the Milky 
Way (MW) seems to exhibit a Holmberg effect even on 
small scales, in that the 11 innermost MW satellites 
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Numerical simulations have been employed to test 
alignment on a similar range in scales, from super- 
clusters down to galaxy-satellite systems. All stud- 
ies focusing on cluster size halos report a corre- 
lation in the orientations for distances of at least 
10/i~^Mpc; some studies observe a positive align- 
_ment signal u p to 100/t ^^Mpc (e.g.. Onuora fc T homas! 
2000 ; Faltcnbacher et al.. 2002, 2005; H opkins et a l.ll200£ 



Kasun fc EvrardI l2005t iBasilakos et all l2006f ). These 
findings are interpreted as the signature of the filamen- 
tary network which interconnects the clusters. The pref- 
erential accretion along these filaments causes the clus- 
ters to point towards each other. Also, for galaxy and 
group-sized halos a tendency to poi nt toward ne i ghbor - 
ing halos is detected. According to lAltay et"an (|2006l ) 
the alignments for such intermediate mass objects are 
caused by tidal fields rather than accretion along the 
filaments. Consequently, the mechanisms responsible 
for the alignment of the orientations depend on halo 
mass. Further evidence for a mass dependence of align- 
ment effects comes fro m the examination of th e ha- 
los' angular momenta. Baffin fc Steinmeti ()2005f ) and 
IXragon-Calvo ct al. (2007) find that the spins of galaxy 
size halos tend to be parallel to the filaments whereas 
the spins of group-sized halos tend to be perpendicular. 
This behavior may originate in the relative sizes of halos 
with respect to the surrounding filaments. 

On subhalo scales basically three different types of 
alignments have been discussed: the alignment of the 
overall subhalo dist ribution with the orientation of 
the host halo (e.g.. iKnebe et al.l l2004t IZentner et al.l 



2005*; 'Agustsson fc BrainerdI I2006U iKang et al.l l2007t 
Libeskind ct al. 2003), the alignme nt of the orienta tions 
of subhalos among each other (e.g.. lLee et al.l l2005') and, 
very recently, the orientation of th e satellites w i th re- 
spect to the cent er of the host (jKuhlen et al.l 120071 : 
iPereira et al.ll2007l ). Again, accretion along the filaments 
and the impact of tidal fields have been invoked as expla- 
nations for the former and the latter, respectively. Thus, 
on all scales tidal fields and accretion along filaments are 
considered to be the main contributers to the observed 
alignment signals. Here we attempt to isolate the dif- 
ferent contributions. In particular we focus on the con- 
tinuous transition from subhalo to halo scales meaning 
we examine the alignment of (sub) structure on distance 
scales between 0.3 and 6 times the virial of groups sized 



Faltcnbacher et al. (2007, hereafter Paper I) 
applie d the halo-based group finder of lYang et all 
(120051) to the SPSS Data Release Four (DR4; 
lAdelman-McCarthy et al"]l200 6^ and carried out a study 
of the mutual alignments between central galaxies (BCG) 
and their satellites in group-sized halos. Using the same 
data set consisting of over 60000 galaxies three different 
types of alignment have been investigated : (1) the 'halo' 
alignment between the orientations of the BCG and as- 
sociated satellite distribution; (2) the 'radial' alignment 
between the direction given by the BCG-satellite con- 
nection line and the satellite orientation; (3) the 'direct' 
alignment between the orientations of the BCG and the 
satellites. The study presented in this paper focuses on 
the same types of alignment and is aimed to compare 
the observational results with theoretical expectations 
derived from N-body simulations. 

There are a variety of dynamical processes which can 
contribute to the alignments of satellites associated with 
groups, the most important are: (1) a possible pre- 
adjustment of satellites in the filaments, which for dis- 
tances of a few times the virial radius commonly point 
radially towards the group; (2) the preferential accre- 
tion along those filaments; (3) the change of the satellite 
orbits in the triaxial group potential well; (4) the contin- 
uous re-adjustment of satellite orientations as they orbit 
within the group. Basically, the first two points can be 
attributed to the large scale environment of the groups 
whereas the latter two are more closely associated with 
the impact of the group potential on small scales. The 
purpose of the present analysis is to separate the dif- 
ferent contributions to the observed alignment signals, 
therefore we analyse the mutual orientations of satellites 
within 6 times the virial radius of the groups. Since the 
tidal forces are closely related to the dynamics of the 
satellites additional insight into the generation of align- 
ment can be gained by considering the satellite velocities. 
Therefore, we also investigate the direction of the satel- 
lite velocities with respect to their orientations, which 
constitutes an indirect way to infer the impact of the 
dynamics onto the orientation of the satellites. A more 
direct way to work out the interplay between the dy- 
namics and the orientations would be to trace the orbits 
of individual satellites, however such an approach goes 
beyond the scope of the present study. 

The paper is organized as follows. In § [2] we introduce 
the simulation and describe the halo finding procedure. 
§[3] deals with some technical aspects, namely the deter- 
mination of the size and orientation of the substructures. 
In § |4] we present the signals of the three dimensional 
spatial and velocity alignments and in § [S] we repeat the 
analysis based on projected data. Finally, we conclude 
with a summary in § [6l 

2. SIMULATION AND HALO IDENTIFICATION 

For the present analysis we employ an N-body sim- 
ulation of structure formation in a flat ACDM uni- 
verse with a matter density ilm = 0.3, a Hubble pa- 
rameter h = i?o/(100 kms~^ Mpc~^) = 0.7, and a 
Harrison- Zeldovich initial power spectrum with normal- 
ization CTg = 0.9. The density field is sampled by 512"^ 
particles within a lOOft-^^Mpc cube resulting in a mass 
resolution of 6.2 x 10^/i~^Mq. The softening length was 
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set to e = lO/i^^kpc, beyond which the gravitational 
force between two particles is exactly Newtonian. The 
density filed is evolved with 5000 time steps from an 
initial redshift of Zi — 72 using a PPPM method. An 
extensive de scription of the simulation can be found in 
iJing fc Sutol p002 ) where it is quoted as LCDMa real- 
ization. 

As detailed in the following two paragraphs the host 
halos and its satellites are found in two subsequent steps 
with two different techniques, first the main halos are 
located thereafter the associated satellite halos are de- 
tected. In orde r to identify the h ost halos we first run a 
FoF algorithm (|Davis et al.lll985l ) on the simulation out- 
put at z = 0. We set the FoF linking length to 0.1 times 
the mean particle separation, which selects regions with 
an average overdensity of ~ 3000. Note that, this link- 
ing length is a factor of two smaller than the commonly 
used value of 0.2, consequently only the central part of 
the host halo (and occasionally large substructures) are 
selected. Subsequently, the virial radius, rvir, is defined 
as the radius of the sphere centered on the most bound 
FoF particle which includes a mean density of 101 times 
the critical density, and we simply define the virial mass 
of each halo as the mass within rvir. If the virial regions 
of two halos overlap, the lower mass halo is discarded. In 
what follows we only focus on the 515 halos with a virial 
mass in the range from 10^^ /i'^Mq to 5 x lO" /i"^Mq 
(corresponding to halos with more than 16,000 particles). 
Since this is the typical mass scale of galaxy groups, we 
will refer to these halos as 'groups'. 

In a second step we search for self -bound 
(sub) structures using the SKID halo finder (jStadell 
1200 It ) applied to the particle distribution within group 
centri c distances of 6 rvir- As discussed in iMaccio et al.l 
(|2006f ) SKID adequately identifies the smallest resolvable 
substructures when using a linking length I equal to 
twice the softening length, i.e. four times the spline 
softening length. We therefore adopt I = 20ft.~^kpc. 
Throughout we will distinguish between "group central 
substructures" (GCSs), which are located at the center 
of our groups, and satellites which are all the other 
(sub)structures, no matter whether they lie within or 
beyond rvir- According to this definition every group 
hosts one, and only one, GCS at its center while it may 
have numerous satellites outside the volume occupied 
by the GCS. Satellites are allocated to all groups from 
which they are separated less than 6 rvir. Hence, a 
satellite may be assigned to more than one GCS. 

3. SIZE AND ORIENTATION OF SUBSTRUCTURES 

Before describing the computation of the orientation 
we determine the typical sizes of the GCSs and the 
satellites. Knowledge about the physical sizes of the 
(sub)structures provides a crucial link for the compar- 
ison to observational data. 

3.1. Sizes of group central substructures 

The physical interpretation of the size of the GCS is 
not straightforward. For one thing, it depends on the 
SKID linking length used. However, for our purposes it 
is sufficient to notice that the GCS represents the dense 
inner region of the group which, largely due to numerical 
reasons, is free of substructure. Consequently, any ra- 
dial dependence of satellite properties can only be probed 



Fig. 1. — Distribution of the spatial dispersion of the group 
central substructure (GCS) in units of the virial radius. Satellites 
can only be resolved at radii larger than the size of the GCS. 
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Fig. 2. — Distribution of the radii of satellites found within 
the virial radius of the group. In this context radius refers to 
listed fraction (0.5, 1.0) of the satellites spatial dispersion ro-. For 
example, the typical inner radii probed by the b = O.Sro- sample is 
~ 30/i~^kpc. 

down to the size of the GCS. In order to express the sizes 
of the GCS and the satellites we use the rms of the dis- 
tances between the bound particles, r^- The advantage 
of this size measure is that it provides a direct estimate 
of the (momentary) size without having to make any as- 
sumption regarding the actual density distribution. In 
the case of an isolated NFW halo r„ « 0.5 rvir, with only 
a very weak dependence on the concentration parameter. 
Figure [1] displays the ra distribution of the GCSs in units 
of the group's virial radius, rvir. The distribution peaks 
at 0.11 rvir and has a mean of 0.13 rvir. 

3.2. Sizes of satellite halos 

The aim of the present analysis is twofold: (1) to as- 
sess the impact of the group tidal field on the satellite 
orientations, and (2) to compare the alignment signals in 
our N-body simulation to observations of galaxy align- 
ments. The impact of the group tidal field is stronger 
at larger satellite-centric radii. On the other hand, since 
galaxies reside at the centers of their dark matter ha- 
los, the central parts of the satellites are more of in- 
terest when comparing the alignment signals with those 
observed for galaxies. To meet both requirements we 
therefore measure the orientation of the satellite mass 
distribution within two radii. In analogy to the mea- 
surement of GCS sizes, we determine these radii with 
reference to the spatial dispersion r„ . More precisely, we 
choose the particles within 1.0 and 0.5rCT as the basic sets 
for the subsequent determination of the satellite orienta- 
tion (see Section [3731 below). Figure [2] displays the dis- 
tributions of the corresponding physical sizes. The 0.5rCT 
sample probes the matter distribution of the satellites 
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Alignments between central and satellite halos 
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Fig. 3. — Distribution of satellite shapes, represented by the ra- 
tio of shortest to the longest semi-major axis of the mass-ellipsoid 
(s = c/a). The colors correspond to the samples with different 
truncation radii as listed. With increasing size the halos become 
rounder. The distribution of the b = O.Sr^ sample is rather sym- 
metric, whereas for larger truncation radii there appears a slight 
asymmetry. 

within ^ 25 /i^^kpc, which is comparable to the sizes of 
elliptical galaxies. The mean, physical radii of the LOr^ 
sample is ~ 50 h~^]ipc. If not quoted otherwise we will 
display the results for the 0.5r„ sample, since this may 
most closely resemble the properties of observable galaxy 
distributions (outside of the very central part of the host 
halo). 

3.3. Orientation 

There are a few diffe r ent w ay s found in th e lit- 
erature (e g.. iBullockl 120021: iJing fc Sutol 
Bailin k Steinmet3 120051 : iKasun fc EvrardI 



2002 



2005 



Allgood et all 120061 ) to model halos as ellipsoids. They 
all differ in details, but most methods model halos using 
the eigenvectors from some form of the inertia tensor. 
The eigenvectors correspond to the direction of the 
major axes, and the eigenvalues to the lengths of the 
semi-m ajor axes a > b > c. Following lAllgood et al.1 
(|200l) we determine the main axes by iteratively com- 
puting the eigenvectors of the distance weighted inertia 
tensor. 

k=l,N I' 



lij — 



(1) 



where r^j denotes the ith component of the position vec- 
tor of the fcth particle with respect to the center of mass 
and 



''fc = A/ — + 52 + ^ 



(2) 



is the elliptical distance in the eigenvector coordinate sys- 
tem from the center to the fcth particle. The square roots 
of the eigenvalues of the inertia tensor determine the ax- 
ial ratios of the halo (a : h : c — ^/\^ : VAb : -s/A^). 
The iteration is initialized by computing the eigenvalues 
of the inertia tensor for the spherically truncated halo. 
In the following iterations the length of the intermediate 
axis is kept unchanged and all bound particles within the 
ellipsoidal window determined by the eigenvalues of the 
foregoing iteration are used for the computation of the 
new inertia tensor. The iteration is completed when the 
eigenvectors have converged. The direction of the result- 
ing major axis is identified as the orientation. The ad- 
vantage of keeping the intermediate axis fixed is that the 
number of particles within the varying ellipsoidal win- 
dows remains almost constant. Instead, if the longest 



(shortest) axis is kept constant the number of particles 
within the ellipsoidal windows can decrease (increase) 
substantially during the iteration. 

Note that we apply this truncation to all 
(sub)structures, both satellites and GCSs, and that 
the orientation of each sub(structure) is measured 
within this truncation radius. Throughout we only 
consider those sub(structures) that comprise at least 
200 bound particles within the volume of the final 
ellipsoid (corresponding to a lower limit in mass of 
« 10^^ /i~^Mq). For the satellites this implies that a 
smaller truncation radius results in a smaller sample. 
For example, there are 772 0.5ro- satellites within the 
virial radii of our groups whereas the LOr^ sample 
comprises 1431 satellites. Since all 515 GCSs contain 
more than 200 particles within 0.5rcr their sample size is 
independent of the truncation radius used. 

Figure [3] displays the distribution of the shape param- 
eter s = c/a. The shading corresponds to different trun- 
cation radii as listed. There is a weak indication that 
satellites become more spherical with increasing trunca- 
tion radi i. A similar behav i or was found for iso lated ha- 
los (e.g.. IJing fc Sutdl2002l:IAllgood et al.ll2f)0i l. As dis- 
cussed bv A llgood et al.l (|2006D the exact determination 
of individual shapes may need as many as 7000 particles, 
so that the resolution of the present simulation is not 
suited for the analysis of (sub)structure shapes. How- 
ever, for the determination of the orientations, which is 
the focus of this paper, a particle l imit of 200 can be con - 
sidered conservative (cf.. lJing|[2002t iPereira et al.ll2007| ). 
A study examining the shapes of substructure in a sin- 
gle high-resol u tion M ilky Way-sized halo can be found in 
iKuhlen et all ([200I . 

4. THREE DIMENSIONAL ALIGNMENTS 

For both classes of objects, GCSs and satellites, the ori- 
entations are determined according to the approach de- 
scribed above. A third orientation-like quantity is given 
by the direction of the line connecting a GCS-satellite 
pair. Throughout we refer to the orientation of the CCS, 
the satellite and the connecting line as accSi ^sat and r, 
respectively. These quantities are unit vectors, such that 
the scalar product of two vectors yields the cosine of the 
angle between them. We will focus on three different 
types of alignment, (1) the halo alignment between the 
orientations of the GCSs and the connecting lines, (2) the 
radial alignment between the orientations of the satellites 
and the connecting lines and (3) the direct alignment be- 
tween the orientation of the CCS and that of its satel- 
lites. In addition, we also consider various alignments 
based on the proper velocity, v, of the satellite with re- 
spect to its GCS. In particular, we discuss (4) the radial 
velocity alignment between v and r, (5) the halo velocity 
alignment between aces and v, and finally (vi) the auto 
velocity alignment between the orientations, agAT, and 
velocities, v, of the satellites. Here v is the unit vec- 
tor indicating the direction of the proper velocity of the 
satellite (including the Hubble flow) relative to the host. 
Since all the other quantities also represent unit vectors 
the scalar products yield the cosines of the enclosed an- 
gles. 

4.1. Halo alignment 
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Fig. 4. — Mean values of the cosines of the angles between the 
orientations of the GCS and the connecting lines to the satellites, 
([aQcs'i"l)> 3- function of r/r^ir for the O.Sro- sample. The 
dotted horizontal line indicates the mean values for an isotropic 
distribution. The error bars indicate the 95% bootstrap confidence 
intervals within each distance bin. 

In order to measure the alignment between the GCS 
and the satellite distribution we use aces ^^nd r (the ori- 
entation of the GCS and the position of the satellite with 
respect to its GCS). Figure H] displays the radial depen- 
dence of (|aGcs'r|) within 6rvir, where (•) denotes the 
mean value within a bin of r/rvir- The error bars in- 
dicate the 95% bootstrap confidence intervals based on 
1000 bootstrap samples for each distance bin. Over the 
entire range of distances probed, the mean values of the 
cosines deviate significantly from a isotropic distribution. 
The strength of the alignment, i.e. the deviation from 
(|aGCS'r|) = 0.5, increases with group centric distance 
and reaches a maximum at ~ l.Trvir- The subsequent 
decline, however, is very weak and even at 6 r^ir the align- 
ment is still very pronounced ((|aGcS'r|) ~ 0.55), with no 
clear indication of a downward trend. The fact that there 
is strong alignment over such a long range suggests that 
the halo intrinsic alignment is closely connected to the 
filamentary structure in which the groups are embedded 
in. Since here we focus on the transition between group 
and environment dominated areas we do not aim to map 
out the entire range of the radial alignment. 

The weakening of the signal at small scales may be 
attributed to the fact that the information about the fil- 
amentary origin is washed away once the satellites start 
to orbit within the groups (i.e., once non-linear effects 
kick in). Yet, the orientation of the group itself is closely 
correlated with the surrounding filamentary network, so 
that a residual alignment is maintained by the overall dis- 
tribution of satellit es orbiting in the p otential well of the 
group (cf. Statler 1 9871 : IZentner et a l. 2005; Kang et al] 
mOTi- Additionally, if one assumes that filaments are ap- 
proximately cylindrical in shape and the GCS is aligned 
with the orientation of the cylinder, then the mean an- 
gles between the orientation of the GCS and the satellites 
position become larger at smaller group-centric radii. In 
fact, at distances smaller than the radius of the cylinder 
the distribution will converge to isotropic. Finally, some 
contribution to the decrease of the alignment strength 
on small scales may come from the fact that satellites on 
nearly radial orbits are filtered out during their epicen- 
ter passage. They get severely stripped and consequently 
the number of particles that remains bound can easily fall 
below the detection criterion (minimum of 200 particles), 
thus weakening the alignment signal. 

4.2. Radial alignment 
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Fig. 5. — Same as Figure|4] but for the distributions of cosines 
between the satellite orientation and the connecting line to the 
GCS, jagAT'i"! for the 0.5 and l.Ora samples. 

The radial alignment, (|agAT'r|) , probes the orienta- 
tions of individual satellites, agAT, relative to the direc- 
tion pointing towards their GCS, r. Figure [5] displays 
(|asAT-r|) for distances up to 6rvir. The line styles rep- 
resent different truncation radii of the satellites. Over 
the entire range of group-centric distances probed, the 
data reveal a significant anisotropic distribution. The 
signal is most pronounced on small scales, where it also 
shows a strong dependence on the truncation radii. The 
l.Orcr sample, which includes the behavior of the outer 
mass shells of the satellites, clearly exhibits a stronger 
deviation from isotropy. Within ~ 1.5 r^ir there is a pro- 
nounced decline of the radial alignment signal, while it 
remains remarkably constant at larger radii. For dis- 
tance in the range between 2 — 6 Tvir we detect a weak 
but significant signal, (laGcS'a-SAxI) ~ 0.52, inconsistent 
with isotropy at 95% co nfidence level in go od agreement 
withlHahn et al.l (|2007D . In a recent study, iKuhlen et al.l 
(Hoof detected no radial alignment for distances ^ 3 Tvir. 
However, their analysis is based on a resimulation of a 
single galaxy-sized host halo. Since this halo is rather iso- 
lated, in that it has not experienced any major merger 
after redshift z — 1.7, it is likely that its ambient fila- 
ments have already been drained. 

At large distances satellites preferentially reside in fil- 
aments (as discussed in the context of Figure S]) which 
point radially towards the groups. Consequently, the sig- 
nal on scales > 2 rvir indicates an alignment between the 
satellite orientations and the filaments in which they are 
embedded. Such an alignment may be caused by accre- 
tion of matter along those filaments or by the local tidal 
fields generated by the mass distribution within the fila- 
ments. The group tidal field is not likely to be responsible 
for the observed large scale alignment signal due to its 
rapid decline with distance. On small scales, however, 
the tidal field can substantia lly alter the orie ntations of 
the satellites. As shown by ICiotti fc Duttal ^9M) the 
time scale on which a prolate satellite can adjust its ori- 
entation to the tidal field of a group is much shorter than 
the Hubble time, but longer than its intrinsic dynamical 
time. Therefore, the adjustment of the satellite orienta- 
tions parallel to the gradients of the group potential offers 
a convincing explanation for a radial alignment signal on 
small scales. This perception is further supported by the 
dependence of the alignment strength on the truncation 
radii of the satellites. For the largest radii, which are 
strongest affected by tidal forces, the alignment signal is 
strongest. 
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Fig. 6. — Same as Figure |4] but for the distributions of cosines 
between the orientation of the satellites and the GCS, laQcs'^sATl- 

4.3. Direct alignment 

The strong signals for halo and radial alignment may 
lead to the expectation of a comparably pronounced sig- 
nal for the direct alignment between the orientation of 
the GCS, aces J and the orientations of its satellites, 
ci-SAT- However, as can be seen in Figure [6l the signal is 
weak. There is only a weak trend for positive alignment 
up to 2rvir. The significance found at distances between 
1 and ~ 2 r^h- seems to be somewhat hi gher 90% con - 
fidence). Based on an analytical model lLee et al.l ()2005f ) 
predict a certain degree of parallel alignment between 
host and satellite orientations due to the evolution of the 
satellites within the tidal shear field of host. The signal 
for the direct alignment may be a relic of this effect. 

To summarize, we find positive alignment signals for 
all three types of alignment tested here. However, they 
differ in strength and radial extent. The halo alignment 
is the strongest and reaches far beyond the virial radii 
of the groups (> 6rvir)- The radial alignment is most 
pronounced at small scales, where it reveals a strong de- 
pendence on the radial extent of the satellite over which 
its orientation has been measured. Although the radial 
alignment is weak beyond ~ 1.5rvir, the signal stays 
remarkably constant out to ~ 6rvir. Finally, the least 
prominent signal comes from the direct alignment. This 
ranking of the alignment strengths is in good agreement 
with the observational results reported in Paper I. 

4.4. Alignments based on subhalo velocities 

If tidal forces give rise to the radial alignment on small 
scales, as displayed in Fig. [5l the satellite orientations 
should be related to their actual velocities and the local 
gradients of the host potential. For instance a satellite 
moving radially towards the GCS will show an enhanced 
radial alignment since the gradient of the potential and 
the actual velocity are pointing in the same direction in- 
ducing an orientation in radial direction. On the other 
side the orientations of satellites moving perpendicular to 
the gravitational field (i.e. tangentially with respect to 
the GCS) will lie in between their velocities and the gra- 
dients of the potential well. To gain some more insight 
into the dynamical origin of the alignments, we include 
the directions of satellite velocities into the alignment 
study. Wc will consider three different kinds of align- 
ments: the radial velocity alignment, |v • r|, the halo 
velocity alignment |aGcs ■ v| and the auto velocity align- 
ment |asAT • v|. To facilitate the interpretation of the 
velocity alignments, we split the satellites according to 
whether their net motion is inward (v • r < 0) or out- 
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Fig. 7. — Ratio of inward moving satellites, /;„ = nin/C^out-l-iin) 
for the O.Sro- sample. Error bars are Poisson. 

ward (v • r > 0) with respect to their group. Figure [7] 
shows the fraction of inward moving satellites, /in, as a 
function of their group centric distances. Note that /jn 
reaches a maximum around ~ 2 rvir , beyond which the 
Hubble flow gradually starts to become more and more 
important. In fact, at sufficiently large radii, where the 
Hubble fiow dominates, one expects that /in = 0, and 
all satellites reveal an outward motion. For satellites 
that are in virial equilibrium within the group potential 
(i.e., at r < Tvir), one expects roughly equal numbers of 
inward and outward moving systems (i.e., /in = 0.5). 
However, on these small scales one has an additional 
contribution from the infall region around the group, 
causing /in > 0.5. In addition, a substantial fraction 
of satellites get stripped below the detection limit (200 
particles) at their peri-centric passage, so that they no 
longer contrib ute to the signal on th eir outward motion 
(cf., Faltcnbac her k Mathews! l2007i ) . At rvir, the out- 
going satellite fraction is about 40%, which is (within 
th e errors) con s istent with the value ~ 30% determined 
bv lWang et all (l2005h . If one assumes an average ratio 
of 6:1 between apo- and peri -center dista nces for typical 
satell ite orbits ((Ghigna et al. 1998; va n den Bosch et al.l 
Il999( l the majority of these satellites m ust have passed 
the c entral parts of the group before fcf.. lDiemand et ahl 
[200l . 

The upper panel of Figure [8] displays the radial velocity 
ahgnmcnt, (|v-r|) , as a function of r/ rvir. (|v-r|) > 0.5 
indicates that the distribution of angles between r and v 
is not isotropic, instead, on average they preferentially 
point in radial directions. This behavior is in agree- 
ment with earlier studies of the velocity anisotropy of 
subhalos which is usually expresse d by the anisotropy 
param eter — l — 0.5{at/ar)^ (e.g.. lBinnev &: Tremaind 
Il987t l. where at and ar denote the velocity dispersions 
of the satellites in the tangential and radial direction, re- 
spectively. Note, (|v-r|) is closely related to f3. If one 
assumes a relaxed (steady-state) halo the above men- 
tioned tendency towards radial motions translates into 
a higher radial velocity dispersion compared to the tan- 
gential one ar > a-fi — at/\/2 (where an and at are the 
one and two dimensional tangential velocity dispersions, 
respectively and tangential isotropy is assumed). Thus 
(|v-r|) > 0.5 on small scales (r < 2 rvir) suggest that 
Cr > at/\/2, in good qualitative agreement with numer- 
ical simulations which have shown that P > for subha- 
los within the virial radius of their hosts. ([Ghigna et al.l 
fl998t iColm eraH[2000t iDiemand et al1[200ll 

In accordance with the spherical collapse model the 
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Fig. 8. — Same as Figure |4] but for (|v-r|) . The upper, middle 
and lower panel displays the signal for all, the inward and the 
outward moving satellites, respectively. 

signal extends out to ~ Zrvir, which roughly reflects the 
distance of turnaround. At 2.5rvir the distribution is 
close to isotropic suggesting that at these distances the 
impact of the group potential is negligible and the satel- 
lite motions are dominated by local potential variations 
arising from the filaments and dark matter halos within 
these filaments. Note that the presence of this filamen- 
tary structure in the vicinity of groups is clearly evident 
from Figure U) Finally, the increase of the radial veloc- 
ity alignment on large scales, > 4rvir, is simply due to 
the Hubble flow (i.e., (|v-r|) ^ 1 at r oo). The 
middle panel of Figure [8] shows (|v-r|) for the inward 
moving satellites only. The radial trend within 2 rvir is 
somewhat enhanced compared to the upper panel. At 
larger radii, the inward moving satellites have a veloc- 
ity structure that is consistent with isotropy. The lower 
panel of Figure [8] reveals a marked difference in the be- 
havior of (|v-r|) for the outward moving satellites. It 
indicates a slightly radial trend for satellites within 1 rvir 
which is much lower than seen in the upper two pan- 
els. Within 1 — 2rvir it drops below 0.5, indicating a 
preference for tangential velocities. Together with the 
information derived from Figure [7] this suggests that a 
substantial fraction of outward moving satellites located 
at 1 — 2 Tvir currently are close to their apo-center pas- 
sage after having crossed the more central regions of the 
group. Finally, on large scales the outward moving satel- 
lites clearly reveal the Hubble flow. 

Figure [9] displays the radial dependence of (laces -vl) 
which measures the cosines of the angels between the 
satellite velocities and the orientation of the GCS. On 
large scales the radial outward motion caused by the 
Hubble flow exceeds the internal velocities of the satel- 
lites within the filaments. Since the GCS is strongly 
aligned with these filaments over the entire radial range 
shown (cf. Figure m, one has that (laces 'vD > 0.5 
on scales where the Hubble flow becomes important 
(> 4rvir)- The strong alignment signal on small scales 
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Fig. 9. — Same as Figure|4l but for (|aGcS'v|> . 
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Fig. 10. — Same as Figure |4l but for the distributions of cosines 
between the satellites velocities and positions, |v-r|. 

indicates that the satellites tend to move parallel to the 
orientation of t he GC S. According to lTorm cn (1997) and 
lAllgood et al.l ()2006( ) the principal axes of the velocity 
anisotropy tensor are strongly correlated with the prin- 
cipal axes of the satellite distribution. Therefore, the 
alignment found for (|aceS'r|) (Figure 2]) actually im- 
plies an analogous signal for (|aces'v|) . However, in 
contrast to the halo afignment, (|accS'r|) , the velocity 
halo ahgnment, (|accs'v|) , only extends out to '--^ 1 rvir- 
Beyond this radius a substantial fraction of the satellites 
shows relatively large angles between their velocities and 
the orientation of the GCS which is consistent with the 
picture of tangential motions associated with the apo- 
center passage of the satellites, as discussed in the con- 
text of Figure [H 

Finally we consider the auto velocity alignment, (|agAT' 
v|) , which reflects the distribution of the cosines between 
the satellite velocities and their orientations, |v • asAxl- 
Fig [To] displays the variation of (|asAT-v|) with rj rvir. 

The signal for (|asAT'v|) shows a maximum at 0.7 rvir. 
At larger distances it decreases quickly. Beyond 1.5 rvir 
it is roughly in agreement with an isotropic distribution. 
A possible reason for the slight central dip is, that satel- 
lites on their peri-center passages move perpendicular to 
the gradients of the group potential. Figure [SJ however, 
revealed a preferential radial orientation of these satel- 
lites. Thus, during the peri-center passages the angles 
between satellite orientations and velocities can become 
large. The degree of the radial alignment depends on the 
ratio between the internal dynamical time of the satel- 
lite, with which it can adjust its orientation, and the 
duration of the peri-center passage. If the peri-center 
passage occurs too quickly the time may be too short 
for a 'perfect' radial alignment (cf.. .Kuhlen et al][2007,) . 
On large scales (1 — 2 rvir) a similar mechanism may take 
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place. Above we have argued that within this distance 
range a substantial fraction of satellites are close to their 
apo-center passage. During this phase the velocities are 
again perpendicular to the gradient of the potential but, 
as indicated by Figure O the satellites are oriented radi- 
ally. The comparison between the signal for (|asAT-r|) 
and (|asAT-v|) suggests that, in a statistical sense, the 
(spatial) radial alignment is maintained during the en- 
tire orbit of the satellite within the potential well of the 
groups, which in turn causes a suppression of (|asAT'v|) , 
at its apo- and pcri-center. 

5. PROJECTED ALIGNMENTS 

To facilitate a comparison with observations, in partic- 
ular with the results presented in Paper I, we repeat the 
foregoing analysis using projected data, i.e. we project 
the particle distribution along one of the coordinate axes 
and compute the second moment of mass for the pro- 
jected particle distribution. Accordingly, for the dis- 
tances between GCS and satellites we use the projected 
values (all satellites within a sphere of 6rvir about the 
GCS are projected), which we label as R (the physical 
distances are labelled as r). 

Since the projections along the three Cartesian coor- 
dinate axes are independent we include all three pro- 
jections of each host-satellite in our 2D sample. To re- 
duce the contamination by satellites associated with mas- 
sive ambient groups we exclude those host-satellite sys- 
tems where another SKID group more massive than the 
GCS (which is most likely the center of an ambient host- 
satellite system) is found within a sphere of 6rvir. After 
rejection of 'contaminated' groups we obtain 1034 and 
543 satellites for the l.Oro- and 0.5rcr samples with 3D 
distances to the GCS < r^-„ (for all groups irrespective 
of their environment we found 1431 and 772, see § 13.31 ) 
Furthermore, since (due to technical reasons) we project 
satellites located within a sphere of 6 r^^^ the projected 
volume at large projected distances shrinks substantially. 
Therefore, we analyze the 2D data only for projected dis- 
tances < 3 Tvir which roughly resembles the projection of 
all satellites within a cylinder with a radius 3 rvir and 
length of lOrvir along the 'line of sight'. Thus, in an ap- 
proximate manner, uncertainties in the determination of 
group membership based on redshift measurements are 
accounted for. 

The resolution of the simulation does not permit to 
probe alignment below 0.3rvir. Other authors (using 
semi-analytical techniques, e.g.. iKang et al.ir2007f ) have 
bypassed this problem by introducing so-called orphan 
galaxies, i.e. galaxies which are associated with the 
once most bound particle of a satellite halo which sub- 
sequently has become undetectable due to the stripping 
by tidal forces. Here we do not adopt this technique 
since it does not provide us with information about the 
orientation of a satellite. Both approaches, considering 
only satellite halos with a minimum number of particles 
and the introduction of orphan galaxies, have certain dis- 
advantages. The former does not account for galaxies 
which are hosted by strongly stripped subhalos whereas 
the latter ignores the dynamical differences of galaxies 
and (once most bound) particles. 

The application of a fixed lower particle limit excludes 
satellites from the analysis which still constitute distinct 
objects. In particular satellites which are strongly tidally 




Fig. 11. — Illustration of the three angles d, <f> and which are 
used for halo al ignment, radial al ignment and direct alignment, 
respectively fcf . . IFaltenbacher et al.ii2007i) . 

stripped may fall below the selection criterion even if the 
galaxy, which is assumed to sit at the center, may still be 
observable. Thus, we caution that our satellite sample 
may be somewhat biased toward more recently accreted 
satellites compared to a hypothetical galaxy population. 
This effect appears whenever a fixed lower particle limit 
is imposed. 

In analogy to Paper I we define the angles 9, (p and ^ to 
address the projected halo, radial and direct alignments 
(same definitions as in Sj4]but for the 2D data, see Fig. fTTj) 
and the projected orientations are referred to as position 
angles (PAs). It is not straightforward to derive galaxy 
properties, such as luminosity and color, from the dark 
matter distribution. In particular, if the satellite halo 
hosts a late type galaxy, it is not obvious how to accu- 
rately determine the orientation of the disk (but s ee e.g., 
iKang et aLll2007l and lAgustsson k Brainerdll2007l for at- 
tempts). On the other side, if one focuses on early type 
galaxies the orientation of the central dark matter distri- 
bution is very likely correlated with the orientation of the 
stellar comp onent (see the e vidence from gravitational 
lensing, e.g., iKochaneS l2002f l . The lower particle limit 
for the satellites results in a lower mass of 10^^ H'^Mq 
within 25/i~^kp c. Assuming a dynam ical mass-to-light 
ratio of a few (jCappellari et al.l l2006f l within this ra- 
dius yields a stellar component which roughly resem- 
bles galaxies. Therefore, our findings in the cur- 
rent paper may be best compared with results based 
on bright early-type satellite galaxies. However, as we 
have pointed out in Paper I, our observational results 
were only marginally dependent on the luminosity/mass 
of satellite galaxies. Therefore, a comparison with obser- 
vations based on somewhat fainter satellites is viable as 
well. 

5.1. Halo alignment 

Figure [T^ shows the results obtained for the angle 6 
between the orientation of the GCS and the line con- 
necting the GCS with the satellite. The short horizontal 
line on the left indicates the result for the innermost bin 
if only the satellites within 1 rvir are projected. The sam- 
ple shows {9) < 45° for the entire distance range. The 
error bars give the 95% bootstrap confidence intervals for 
the mean angles within each bin. The alignment strength 
wi thin Tyir i s ^ 42 °, in good agr e ement with the findings 
of iBrainerdI (|2005D . lYang et al.l (|2006[ ). In Paper I we 
found a mean value 9 « 41° within 0.5 rvir which is very 
close to the values we obtain for the innermost bin, in 
particular if only the satellites within 1 rvir (short hori- 
zontal lines on the left) are projected. As also shown by 
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Fig. 12. — Mean angle, 6*, between the PA of the GCS and the 
Une connecting the GCS and a sateUite, as a function of projected 
distance R/ rvir with equidistant bins of 0.5 X rvir- The error bars 
give the 95% bootstrap confidence intervals for the mean angles 
within each bin. The short horizontal line on the left indicates 
the signal for the innermost bin if only the satellites with in 3D 
distances < 1 rvir are projected. The corresponding 3D results are 
shown in Fig. |4] 
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Fig. 13. — Same as Figure [T2l but for the angle </>. In addition 
the radial dependence of the b = Ir^ sample is displayed as well. 
The corresponding 3D results are shown in Fig. |5] 

lAgustsson fc BrainerdI (|2006H ) the alignment signal ex- 
tends beyond the virial radius. The strongest amplitude 
is found outside the virial radius at ~ 1.7 rvir- Currently 
there are no available observations covering the same dis- 
tance range. The analysis in Paper I, for instance, is 
based on galaxies within the virial radius whereas we use 
all galaxies with projected distances < 3 rvir- According 
to our findings a search for alignment of satellite distri- 
bution in group environments for distances larger than 
rvir may be a promising proposition. 

5.2. Radial alignment 

Figure [13] displays the mean angle </> between the PA of 
the satellite and the line connecting the satellite with its 
GCS. For all group centric distances there is a clear and 
significant signal for the major axes of the satellites to 
point towards the GCS (i.e., {((>) < 45°). The projection 
of only those satellites within 1 rvir increases the central 
signal by about 1° (differences between the innermost 
data points and the short horizontal lines). The mean 
angle for the O.Sro- sample within the innermost bin is 
~ 43° and according to Paper I the mean value for the 
red SDSS satellites within 0-5 rvir is very close to this 
value. However, the observations suggest a significant 
alignment for red galaxies only out to 0.7 rvir whereas 
the N-body data indicate that radial alignment extends 
beyond 3 rvir- The discrepancy may be caused by the 
observational confinement to galaxies within the virial 
radius. 
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Fig. 14. — Same as Figure [T2l but for the angle 

5.3. Direct alignment 

Figure [13] displays the results for the direct alignment, 
based on the angle ^ between the orientations of GCSs 
and satellites. The alignment signal is significant at a 
> 95% confidence level for distances < 0.5 rvir- In Pa- 
per I we obtained ^ « 44° for red satellite with in 0.5 rvir 
which indicates a somewhat weaker alignment than we 
find here. Since the 3D analysis shows no increase of 
(|aGcs'r|) at small scales (Figure HJ the central enhance- 
ment displayed here has to be interpreted as a result of 
projection effects. 

In summary for all three types of alignments we find 
good agreement between numerical data presented here 
and the observational results from Paper I. In particu- 
lar the relative strength among the different alignments is 
well reproduced in the numerical analysis. Due to limited 
resolution the range below 1 rvir is only sparsely sampled 
thus no detailed information about the radial dependence 
of the alignment signal on small scales can be derived. 
However, the signal for increases with distance which 
is only marginally implied by the SDSS results presented 
in Paper I. Also for (j), the dependence on the distance 
disagrees between simulations and observations. It is cur- 
rently unclear whether this is due to shortcomings from 
the numerical or observational side. 

6. SUMMARY 

Based on a sample of 515 groups with masses rang- 
ing from 10^^ /i^^Mq to 5 x 10^*^ H'^Mq we have inves- 
tigated the halo alignment, (|aGcs'r|) , the radial align- 
ment, (|asAT'r|) and the direct alignment (laGcS'^sATl) , 
between the central region of each group (the GCS) and 
its satellite halos out to a distance of 6 rvir- Here aces, 
a-SAT and r denote the unit vectors associated with the 
orientation of the GCS, the satellites and the line con- 
necting both of them. Additionally, we have employed 
the directions of the satellite velocities v to probe the 
ahgnments (|v-r|) , (|aGcS'v|) and (|asAT-v|), referred 
to as radial, halo and auto velocity alignments, respec- 
tively. Our main results are: 

(1) Halo, radial and direct alignment differ in strength. 
The halo alignment is strongest followed by the ra- 
dial alignment. By far the weakest and least signif- 
icant signal comes from the direct alignment. This 
sequence is found in the 3D analysis as well as for 
the projected data and agrees well with our recent 
analysis of galaxy alignments in the SDSS (cf.. Pa- 
per I). 

(2) The signal for the halo alignment, (|aGcs ■ r|) , 
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reaches far beyond the virial radii of the groups 
(> 6rvir) which we interpret as evidence for large 
scale filamentary structure. 

(3) The signal for the radial alignment, (|agAT'r|) , is 
largest on small scales. After a rapid decline with 
distance it flattens, such that a relatively small 
(|asAT-i"|) ~ 0.52, but significant deviation from 
isotropy is detected out to ~ 6rvir. Whereas the 
small scale signal more likely owes to the group's 
tidal field, the weak but significant signal on large 
scales suggests that satellites tend to be oriented 
along the filaments in which they reside. 

(4) The 3D signal for the direct alignment, (|aGcs • 
agAxI) : shows a weak trend for parallel orientations 
on scales < 2rvir. The projected data indicate an 
increasing signal for distances < 0.5 rvir which is 
likely caused by projection effects. 

(5) All kinetic alignment signals are highly significant 
at small scales. The signal for (|v-r|) is basically 
constant within 2.0 rvir, beyond which it rapidly 
drops. In the subset of outward moving satellites 
we find a tendency for tangential motions which 
can be attributed to the satellites which have been 
accreted earlier and are currently passing their 
peri- or apo-centers. The signal for (|aGcS'v|) is 
maximal at the center, drops rapidly with distance 
and disappears at Irvir- Finally, (|asAT'v|) shows 
a slight dip at the center, reaches a maximum at 
0.7rvir, and becomes consistent with isotropy at 
1.5 rvir- AH these features support the interpreta- 
tions advocated for the spatial alignments. 

The simulation analyzed here clearly demonstrates that 
tidal forces cause a variety of alignments among neigh- 
boring, non-linear structures. On large scales, the tidal 
forces are responsible for creating a filamentary network, 
which gives rise to a halo alignment out to at least 6 rvir- 
The same tidal forces also cause an alignment between 
filaments and (sub)structures within the filaments (cf . , 
lAltav et all l2006l : iHahn et all l2007t l which in turn re- 
sults in a large scale radial alignment with the virialized 
structures at the nodes of the cosmic web. Within these 
virialized structures, tidal forces are responsible for a ra- 
dial alignment of its substructures, similar to the tidal 
locking mechanism that affects the Earth-Moon system. 
This is further supported by the fact that the auto ve- 
locity ahgnment (|agAT-v|) reveals a dip on small scale, 
indicating that at peri-centric passage satellites tend to 
be oriente d perpendicu l ar t o the direction of their mo- 
tion (cf., Kuh len et al.1 '2007). This behavior also ex- 
plains, why the direct spatial alignment, (|aGcs 'a-SATl) , 
is so weak. A possible direct alignment originating from 
the co-evolution of group and satellites, as proposed by 
iLee et aP ()2005| ). is quickly erased as the satellites orbit 
in the potential well of the group. For future work it will 



be instructive to trace the orbits of individual satellites 
and consider more closely how their shapes and orienta- 
tions evolve with time. 

The infall regions around virialized dark matter ha- 
los cause a radial velocity alignment out to ^ 2 rvir, 
and an enhancement of inward moving (sub)structures. 
At around the same scale, the (sub)structures with a 
net outward movement have a tendency to move tan- 
gentially. This most likely reflects the apo-centric pas- 
sage of substructures that have previously fallen through 
the virialized halo. Within a virialized region, the ori- 
entation of orbits is naturally aligned with that of its 
GCS. Since (sub)structures reveal at most a weak ve- 
locity bias with respect to da rk matter particles (e.g., 
iFaltenbac her fc Diem and 20061). this causes a strong halo 
velocity alignment on scales < rvir. The halo veloc- 
ity alignment is also strong on large scales (> 3 rvir), 
which reflects the Hubble flow combined with the fila- 
mentary, non-isotropic distribution of (sub)structures on 
these scales. 

A one-to-one comparison between the N-body results 
discussed here and the observations presented in Pa- 
per I is not straightforward. Although we have employed 
the same mass range for the groups in both studies the 
resolution of the current simulation only allows to re- 
solve satellites which are expected to host > galaxies. 
These are bright compared to our SDSS sample for which 
a lower magnitude limit of ^-^Mr — 5 log ft. < —19 has 
been adopted. Nevertheless, the qualitative agreement 
between the relative strengths of the different types of 
spatial alignment is promising. Supplementary to the 
observational results of Paper I we find a strong halo 
alignment and a somewhat weaker radial alignment out 
to at least 6 rvir which we will investigate further. 

Finally, the weak but significant detection of radial 
alignment out to 6 rvir may contaminate the cosmic shear 
measurements on these scales. This correlation has 
to be considered, either by simply removing or down- 
weighting pairs of galax ies within this distance range 
(|King fc Schneided l2002l : iHevmans fc Heavens! l2003f i. 
This may be particularly important for applications of 
weak gravitational lensing for the purposes of precision 
cosmology. 
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